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^ I An efficient technique to generate ensembles of spins that are 

O I highly polarized by external magnetic fields is the Holy Grail in 

Nuclear Magnetic Resonance (NMR) spectroscopy. Since spin-half 
nuclei have steady-state polarization biases that increase inversely 
! with temperature, spins exhibiting high polarization biases are con- 

I sidered cool, even when their environment is warm. Existing spin- 

cooling techniques [V, 4j are highly limited in their efficiency 

and usefulness. Algorithmic cooling [5\ is a promising new spin- 
, , cooling approach that employs data compression methods in open 

'jjj . systems. It reduces the entropy of spins on long molecules to a 

I point far beyond Shannon's bound on reversible entropy manipula- 

tions (an information-theoretic version of the 2nd Law of Thermo- 
dynamics), thus increasing their polarization. Here we present an 
efficient and experimentally feasible algorithmic cooling technique 
that cools spins to very low temperatures even on short molecules. 
This practicable algorithmic cooling could lead to breakthroughs 
in high-sensitivity NMR spectroscopy in the near future, and to 
the development of scalable NMR quantum computers in the far 
future. Moreover, while the cooling algorithm itself is classical, it 
uses quantum gates in its implementation, thus representing the 
first short-term application of quantum computing devices. 
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NMR is a technique used for studying nuclear spins in magnetic fields jS]. 
NMR techniques are extremely useful in identifying and characterizing chem- 
ical materials, potentially including materials that appear at negligible lev- 
els. There are numerous NMR applications in biology, medicine, chemistry, 
and physics, for instance magnetic resonance imaging (used for identifying 
malfunctions of various body organs, for monitoring brain activities, etc.), 
identifying materials (such as explosives or narcotics) for security purposes, 
monitoring the purity of materials, and much more. A major challenge in 
the application of NMR techniques is to overcome difficulties related to the 
Signal-to-Noise Ratio (SNR). Several potential methods have been proposed 
for improving the SNR, but each of them has its problems and limitations. 
We thoroughly review in the Supplementary Information A the six main 
existing solutions to the SNR problem. In brief, the first (and not very effec- 
tive) three methods involve cooling the entire system, increasing the magnetic 
field, and increasing the sample size. The fourth method involves repeated 
sampling over time, a very feasible and commonly practised solution to the 
SNR problem. Its limitation is that in order to improve the SNR by a factor 
of M, spectroscopy requires repetitions, making it an overly costly and 
time-consuming method. 

The fifth and sixth methods for improving the SNR provide ways for 
cooling the spins without cooling the environment, which in this case refers 
to the molecules' degrees of freedom. This approach, known as "effective 
cooling" of the spins [H 121 E], is at the core of this current work, and is 
explained in more detail in Supplementary Information A. The spins cooled 
by the effective cooling can be used for spectroscopy as long as they have 
not relaxed back to their thermal equilibrium state. For two-level systems 
the connection between the temperature, the entropy, and the population 
probability is a simple one. The population difference between these two 
levels is known as the polarization bias. Consider a single spin-half particle 
in a constant magnetic field. At equilibrium with a thermal heat bath the 
probability of this spin to be up or down (i.e. parallel or anti-parallel to the 
field's direction) is given by: P| = and P| = The polarization bias, 
e = P| — P|, is around 10~^ — 10~^ in conventional NMR systems, so all the 
following calculations are done to leading order in e. A spin temperature at 
equilibrium is T = ^^S2S^_ \ single spin can be viewed as a single bit with 
"0" meaning spin-up and "1" meaning spin-down, and thus, the single spin 
(Shannon's) entropy is H = 1 — (e^/ In 4). A spin temperature out of thermal 
equilibrium is still defined via the same formulas. Therefore, when a system 
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is removed from thermal equilibrium, increasing the spins' polarization bias is 
equivalent to cooling the spins (without cooling the system) and to decreasing 
their entropy. 

One of these latter two methods for effective cooling of the spins is the 
reversible polarization compression (RFC). It is based on entropy manipula- 
tion techniques (very similar to data compression), and can be used to cool 
some spins (bits) while heating others j2[ E] ■ Contrary to conventional data 
compression, RPC techniques focus on the low-entropy bits (spins), namely 
those that get colder during the entropy manipulation process. The other 
effective cooling method is known as polarization transfer Ij: If at thermal 
equilibrium at a given temperature, the spins we want to use for spectroscopy 
(namely, the observed spins) are less polarized than other spins (namely, aux- 
iliary spins), then transfering polarization from the auxiliary highly polarized 
spins into the observed spins is equivalent to cooling the observed spins, while 
heating the auxiliary spins. In its most general form the RPC can be ap- 
plied to spins which have different initial polarization biases. It follows that 
polarization transfer is merely a special case of RPC. Therefore, we refer to 
these two techniques together as RPC. 

RPC can be understood as reversible, in-place, lossless, adiabatic entropy 
manipulations in a closed system. Unfortunately, as in data compression, 
RPC techniques are limited in their cooling applicability because Shannon's 
bound states that entropy in the closed system cannot decrease. [Shannon's 
bound on entropy manipulations is an information-theoretic version of the 
2nd Law of Thermodynamics.] Polarization transfer is limited because the 
increase of the polarization bias is bounded by the spin polarization bias of 
the auxiliary highly-polarized spins. RPC done on n uncorrelated spins with 
equal biases requires extremely long molecules in order to provide significant 
cooling. The total entropy of such spins satisfies H{n) = n{l — e^/ln4). 
Therefore, with e = 10~^ for instance, molecules of length of an order of 
magnitude of 10^° are required in order to cool a single spin close to zero 
temperature. For more modest (but still significant) cooling, one can use 
smaller molecules (with n <^ 10^°) and compress the entropy into n — 1 fully- 
random spins. The entropy of the remaining single spin satisfies if (single) > 
1 — ne^/ In 4, thus we can at most improve its polarization to 



Figure CJA. provides an illustration of how Shannon's bound limits RPC tech- 
niques. Unfortunately, manipulating many spins, say n > 100, is a very 
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difficult task, and the gain of -y/n in polarization here is not nearly substan- 
tial enough to justify putting this technique into practice. 

We and our colleagues (Boykin, Mor, Roychowdhury, Vatan and Vrijen, 
hereinafter refered to as BMRVV), invented a promising new spin-cooling 
technique which we call Algorithmic cooling 0. Algorithmic cooling ex- 
pands the effective cooling techniques much further by exploiting entropy 
manipulations in open systems. It combines RPC with relaxation (namely, 
thermalization) of the hotter spins, as a way of pumping entropy outside 
the system and cooling the system much beyond Shannon's bound. In order 
to pump the entropy out of the system, algorithmic cooling employs regu- 
lar spins (which we call computation spins) together with rapidly relaxing 
spins. Rapidly relaxing spins are auxiliary spins that return to their ther- 
mal equilibrium state very rapidly. We refer to these spins as reset spins, or, 
equivalently, reset bits. The ratio, -Rreiax-times, between the relaxation time of 
the computation spins and the relaxation time of the reset spins must satisfy 
-Rreiax-times ^ 1- Tliis is vital if wc wish to perform many cooling steps on 
the system. 

The BMRVV algorithmic cooling used blocks of m computation bits and 
pushed cooled spins to one end of such blocks in the molecule. To obtain 
significant cooling, the algorithm required very long molecules of hundreds 
or even thousands of spins, because its analysis relied on the law of large 
numbers. As a result, although much better than any RPC technique, the 
BMRVV algorithmic cooling was still far from having any practical implica- 
tions. 

In order to overcome the need for large molecules, and at the same time 
capitalize on the great advantages offered by algorithmic cooling, we searched 
for new types of algorithms that would not necessitate the use of the law of 
large numbers for their analysis. Here we present examples of a novel, ef- 
ficient, and experimentally feasible algorithmic cooling technique which we 
name ^^practicable algorithmic cooling" . The space requirements of practica- 
ble algorithmic cooling are much improved relative to the RPC, as we illus- 
trate in FiguresQ^ and^C. In contrast to the BMRVV algorithm, practicable 
algorithmic cooling techniques determine in advance the polarization bias of 
each cooled spin at each step of the algorithm. We therefore do not need to 
make use of the law of large numbers in the analysis of these techniques, thus 
bypassing the shortcomings of the BMRVV algorithmic cooling. Practicable 
algorithmic cooling allows the use of very short molecules to achieve the same 
level of cooling as the BMRVV achieves using much longer molecules, e.g.. 
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34 spins instead of 180 spins, as we demonstrate explicitly in Supplementary 
Information B. 

Both RFC and algorithmic cooling can be understood as applying a set 
of logical gates, such as NOT, SWAP etc., onto the bits. For instance, one 
simple way to obtain polarization transfer is by a SWAP gate. In reality, 
spins correspond to quantum bits. Here we provide a simplified "classical" 
explanation, and the more complete "quantum" explanation is provided in 
Supplementary Information C. A molecule with n spins can represent an 
n-bit computing device. A macroscopic number of identical molecules is 
available in a bulk system, and these molecules act as though they are many 
computers performing the same computation in parallel. To perform a de- 
sired computation, the same sequence of external pulses is applied to all the 
molecules/computers. 

Algorithmic cooling is based on combining three very different operations: 

1. RFC steps change the (local) entropy in the system so that some com- 
putation bits are cooled while other computation bits become much 
hotter than the environment. 

2. Controlled interactions allow the hotter computation bits to adiabat- 
ically lose their entropy to a set of reset bits via polarization transfer 
from the reset bits into these specific computation bits. 

3. The reset bits rapidly return to their initial conditions and convey their 
entropy to the environment, while the colder parts (the computation 
bits) remain isolated, so that the entire system is cooled. 

By repeatedly alternating between these operations, and by applying a re- 
cursive algorithm, spin systems can be cooled to very low temperatures. 

In our simplest practicable algorithm there is an array of n computation 
bits, and each computation bit X has a neighbouring reset bit rx with which 
it can be swapped. The algorithm uses polarization transfer steps, reset steps 
and the 3-bit RFC step which we call 3-bit-compression (3B-Comp). The 
3B-Comp can be built, for instance, via the following two gates: 1. — Use 
bit C as a control, and bit 5 as a target; apply a CNOT (ControUed-NOT) 
operation: C — > C, S — > S © C, where © denotes a logical eXclusive OR 
(XOR). 2. — Use bit B as a control, and bits A and C as targets; apply a 
variant of a C-SWAP operation: A ^ CB + AB, B ^ B, C ^ AB + CB; 
this means that A and C are swapped ii B = 0. The effect of a 3B-Comp 
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step in case it is applied onto three bits with equal bias e -C 1 is to increase 
the bias of bit A to 

^new = — • {^) 

See Supplementary Information D for calculation details. A simplified exam- 
ple for practicable algorithmic cooling based on these three steps and using 
77, = 3 computation bits is presented in Figure |21 

In many realistic cases the polarization bias of the reset bits at thermal 
equilibrium is higher than the polarization bias of the computation bits. By 
performing an initiation process of polarization transfer from relevant reset 
bits to relevant computation bits, prior to any polarization compression, we 
cool the computation bits to the initial bias of the reset bits, eo, the zero'th 
purification level. See for instance Figure El 

In the algorithm we now present, we attempt to cool a single bit much 
more rapidly than possible in the reversible processes. In order to cool a 
single bit (say, bit A) to the first purification level, ei = 3eo/2, start with 
three computation bits, ABC, perform a Polarization Transfer (PT) step to 
initiate them, and perform 3B-Comp to increase the polarization of bit A. In 
order to cool one bit (say, bit A) to the second purification level (polarization 
bias €2), start with 5 computation bits {ABODE). Perform a PT step in 
parallel to initiate bits ABC, followed by a 3B-Comp that cools A to a bias 
ei. Then repeat the above (PT + 3B-Comp) on bits BCD to cool bit B, 
and then on bits CDE to cool bit C as well. Finally, apply 3B-Comp to 
bits ABC to purify bit A to the second purification level 62, which for small 
biases gives €2 ~ (3/2)ei ~ (3/2)^eo. 

This (very simple) practicable algorithmic cooling, which uses the 3B- 
Comp, PT, and RESET steps, is now easily generalized to cool a single 
spin to any purification level J/. Let j be an index telling the polarization 
level j G 1 . . . J/. To obtain one bit at a purification level j we define the 
procedure Mj. Let Mq{A) be the PT step from a reset bit to bit A to yield 
a polarization bias eo. Then the procedure Mi contains three such PT steps 
(performed in parallel on three neighbouring bits) followed by one 3B-Comp 
that cools the left bit (among the three) to level 1. In order to keep track 
of the locations of the cooled bit we mention here the location of the cooled 
bit. For an array of n bits, a„a„_i . . . 0201 we write Mi{k) to say that Mi is 
applied on the three bits a^; dk-i] O'k-2 so that bit is cooled to a bias level 1. 
Similarly, M2{k) means that M2 is applied on the five bits a^; ak-i, ■ ■ ■', 0^-4, 
so that bit is cooled to a bias level 2. We use the notation S{(j_i)^j}(fc) 
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to present the 3B-Comp applied onto bits a^] ak-i] ak--2i purifying bit 
from ej„i to tj. Then, the full algorithm has a simple recursive form: for 

je{l,...,J/} 

M,{k) = i3{(,_i)_,}(fc)M,_i(fc - 2) M,_i(fc - 1) M,_i(A;) , (3) 

applied from right to left (Mj_i(fc) is applied first). For instance, Mi (3) = 
S{o^i}(3)Mo(l) Mo(2) Mo(3) is the 3B-Comp applied after reset. The pro- 
cedure cooling one bit to the second level (starting with five bits) is written 
as M2(5) = i3{i^2}(5)Mi(3) Mi (4) Mi (5). Clearly, Mi(fc) can be applied to 
any k >3, M2{k) can be applied to any > 5, and Mj{k) can be applied to 
any k>2j + l. Thus, if we wish to bring a single bit to a polarization level 
J/, then 2Jf + 1 computation bits and 2J/ + 1 reset bits are required. We 
conclude that the algorithm uses 

AJf + 2 (4) 

bits. Due to Eq.([2I) the final polarization bias is egnai = ^ J/ ~ (3/2)'^^eo, as 
long as ej^ -C 1. 

Let us now calculate the time complexity of the algorithm. The number 
of time steps in operation Mj presented above is Tj with Tq = 1, Ti = 2 
(since the three resets are done in parallel), T2 = 7, and Tj = 3Tj_i + 1, for 
any j > 2. The recursive formula yields Tj = 3^~^Ti + (3-'^^ — l)/2, so that 
finally 

5 X - 1 

y 

If we count the reset steps only, assuming a much longer time for reset steps 
than for the other steps, then we get Tj^ (reset) = S'^f^ reset time steps. In 
Supplementary Information B we compare the space and time complexity of 
this practicable algorithmic cooling with the BMRVV algorithmic cooling. 

If we use the reset bits also for the computation, we can improve the space 
complexity by a factor of 2. The algorithm then uses 2Jf computation bits 
and one reset bit, thus a total of 2J/ + 1 bits. See Supplementary Information 
E for a full description of this algorithm. In Figure Q we compare the two 
algorithms descibed here with the RFC, to explicitly illustrate the advantages 
of practicable algorithmic cooling. 

We suggested here an algorithmic cooling technique which is feasible 
with existing technology. A variant of the 3B-Comp step was implemented 
by [H] following Schulman and Vazirani's RFC (their molecular-scale heat 



Tjf = ^ • (5) 



7 



engine) 0. The PT plus RESET steps have been implemented by us and 
our colleagues [3 |H] following the theoretical ideas presented in this cur- 
rent work. Thus, all steps needed for implementing practicable algorithmic 
cooling of spins have already been successfully implemented in NMR lab- 
oratories. Both these implementations use Carbon molecules that contain 
two enriched ^^C atoms which have spin-half nuclei. In addition to the Car- 
bons, the 3B-Comp experiment ^ also uses one Bromine nucleus, the 3-bit 
molecule being C2F3Br. The compression yields a single cooled bit with a 
polarization bias increased by a factor of 1.25. The experimental PT and 
RESET steps j?! |H] use one Hydrogen and the two Carbons, and the 3-bit 
molecule is the Trichloroethylene (TCE) molecule (see Figure EI), with the 
two Carbons being the computation bits and the Hydrogen being the reset 
bit. This work of [71IH] implements our ideas in order to prove experimentally 
that it is possible to bypass Shannon's bound on entropy manipulations. It 
is important to note that cooling spins to high purification levels is a feasible 
task but certainly not an easy one. Addressing specific spins in molecules 
containing many spins is a very challenging operation from an experimental 
point of view. Another important challenge is achieving a sufficiently good 
ratio -Rreiax-times that wiU allow the performance of many reset operations, 
before the cooled bits naturally re-thermalize. 

As demonstrated here the theoretical practicable algorithmic cooling is a 
purely classical algorithm. However, the building blocks are spin-half par- 
ticles which are quantum bits, and NMR machine pulses which implement 
quantum gates, that have no classical analogues. Practicable algorithmic 
cooling is, therefore, the first short-term application of quantum comput- 
ing devices'^ using simple quantum computing devices to improve SNR in 
NMR spectroscopy and imaging. Both practicable algorithmic cooling and 
quantum key distribution could lead to implement 

Algorithmic cooling could also have an important long-term application, 
namely, quantum computing devices that can run important quantum al- 
gorithms such as factorizing large numbers ^U]. NMR quantum comput- 
ers [HI 1121 currently the most successful quantum computing devices 
(see for instance, ^31); but are known to suffer from severe scalability prob- 
lems [TU[T31|3]. Impressive theoretical solutions were provided in [3 E] but 
only the algorithms presented here can lead to a realistic solution. As ex- 

^ Practicable algorithmic cooling is also one of the first short-term applications of quan- 
tum information processing 
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plained in Supplementary Information Al, practicable algorithmic cooling 
can be used for building scalable NMR quantum computers of 20-50 quan- 
tum bits if electron spins will be used for the PT and RESET steps. 
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Figure 1: Comparing space complexity of RFC and practicable algorithmic 
cooling. We present here the number of spins required in order to improve the 
polarization bias of a single spin 5 times via RFC and via practicable algo- 
rithmic cooling 1 and 2, when all spins have the same initial polarization bias 
e. The colors illustrate the temperature, black being the thermal equilibrium 
temperature, blue being colder (five times colder, in the case we present), 
and red being hotter. Case A (RFC): the number of required spins is calcu- 
lated from Shannon's bound (Eq^) — 25 spins are required. [Note that any 
specific algorithm shall probably be much worse than this theoretical bound, 
requiring many more spins.] Case B (FACl): the algorithm "FACl" (for 
Fracticable Algorithmic Cooling 1) is the algorithm described in the text, 
and summarized in Eq. |3] — 18 spins are required. Note that with one addi- 
tional single (parallel) FT step and one additional RESET step all hot spins 
regain their initial bias. Case C (FAC2): the algorithm "FAC2" (for Frac- 
ticable Algorithmic Cooling 2) is described in Supplementary Information E 
— 9 spins are required. [A remark. In order to improve the polarization 25 
times the numbers are: RFC — 625 spins are required; FACl — 34 spins 
are required; FAC2 — 17 spins are required. We did not provide figures due 
to the difficulty of providing the figure for the RFC case.] 
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Figure 2: A molecule for running a simplified algorithm. We show here an 
abstract example of a molecule with three computation bits ABC, such that 
each is attached to a reset bit. All these bits have the same polarization bias. 
The temperature after each step is illustrated by the colors: yellow means 
the temperature at thermal equilibrium; blue means colder than the initial 
temperature; and red means hotter than the initial temperature. Let three 
computation bits A, B, and C be given, such that each is attached to a reset 
bit {ta, rB, and rc), and all bits have the same bias eo (see FigureEK)- Our 
algorithm uses a particular RFC step on 3 bits that we call 3-bit compres- 
sion (3B-Comp), Folarization Transfer (FT) steps, and RESET steps: 1. — 
3B-Comp(yl; S; C); the outcome of this step is shown in Figure Eb- 2. — 
FT(rB B), FT(rc — > C); the outcome of this step is shown in Figure |2fc. 
3. — RESET(rB, rc); the outcome of this step is shown in Figure Eli. The 
3-bit-compression used in the first step operates on the three computation 
bits, increasing the bias of bit A by 3/2 times, while heating the other two 
spins. Note that in the last step also relaxes, but is not modified by the 
RESET operation. As a result of the 3B-Comp step that cools bit A, and 
the later FT and RESET steps that cause all the other bits to regain their 
initial biases, the entire system is cooled down. 
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Tri-Chloro-Ethylene (TCE) 




Figure 3: A 3-bit computer: a TCE molecule labeled with two C. The TCE 
molecule has three spin-half nuclei, the two Carbons and the Hydrogen, which 
are given the names CI, C2 and H, as portrayed in the figure. Only these 
spins take a non-negligible part in the process, so the TCE molecule can 
therefore act as a three-bit computer. The Chlorines have a very small signal 
and their coupling with the Carbons can be averaged out so they cannot be 
considered as bits. TCE can be used for polarization transfer because the 
Hydrogen's nucleus is four times more polarized than the Carbons' nuclei. 
The Hydrogen can also be used as a reset bit, because its relaxation time can 
be made much shorter than the relaxation time of the Carbons. Based on 
the theoretical ideas presented in this current work, we and colleagues [3 |H1 
used the Hydrogen in the TCE to cool both spins CI and C2, decrease the 
total entropy of the molecule, and bypass Shannon's bound on cooling via 
RFC. 
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Supplementary Material 



A Solutions to the SNR problem and their 
limitations 

There are several methods which provide rather trivial solutions to the SNR 
problem. The first of these involves simply cooling the entire system, in- 
cluding cooling the spins, thereby improving the SNR. This is not usually 
a useful solution because cooling the system can result in modifying the in- 
spected molecules (e.g., solidifying the material), or destroying the sample 
(e.g., killing the patient). A second method involves increasing the magnetic 
field, which in some sense is equivalent to cooling the system. This solution 
is generally much too expensive to be considered feasible, as it adds tremen- 
dously to the cost of the necessary machinery. The third method involves 
increasing the sample size, such as taking a larger blood sample. It is usually 
impossible, however, to implement this solution due to the large sample size 
that would be required in order to improve the SNR. One would need to 
increase the sample size by a factor of in order to improve the SNR by 
a factor of M. In the case of a blood sample, then, to obtain a one hun- 
dredfold improvement in SNR one would need to take 100,000 ml. of blood 
instead of 10 ml., which is clearly not feasible. In some cases increasing the 
sample size is strictly impossible, for example if the sample in question is 
a human brain. NMR machine size limitations also severely limit the ap- 
plicability of this method. The fourth method involves repeated sampling 
over time. This is a very feasible and commonly practised solution to the 
SNR problem. Its limitation is that in order to improve the SNR by a factor 
of M, spectroscopy requires repetitions, making it an overly costly and 
time consuming method. In many cases it is also simply impractical, due 
to changes over time in the sample during the spectroscopy process (e.g., 
changes in concentration of a material in a certain body organ). The long 
recovery time between measurements also contributes to the impracticality of 
this method. Furthermore, this solution is less useful if the noise in question 
is not a Gaussian noise. 

An additional two methods provide ways for cooling the spins without 
cooling the environment, which in this case refers to the molecules' degrees 
of freedom. This approach is known as "effective cooling" of the spins. Such 
cooling is as good as regular cooling of the system, because the cooled spins 
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can be used for spectroscopy as long as they have not relaxed back to their 
thermal equilibrium state. For k-Icvcl quantum systems there is, generally, a 
complex relationship between temperature, entropy, and population density 
at the different levels. For the sake of simplicity we shall focus here only 
on two-level systems, namely spin-half particles. Note, however, that the 
results here can be generalized also to higher level systems, and to quantum 
systems other than spins. For two-level systems the connection between the 
temperature, the entropy, and the population density is a simple one. In 
two-level systems the population density difference is known as the polariza- 
tion bias. Cooling the system and therefore also the spins (after they reach 
thermal equilibrium) results in increasing the spins' polarization bias, and 
decreasing their entropy. Consider a single spin in a constant magnetic field. 
At equilibrium with a thermal heat bath the probability of this spin to be 
up or down (i.e. parallel or anti-parallel to the field's direction) is given by: 
P| = ij^, and = The polarization bias, e — — Pi, is calculated 

to be e = tanh (2^^^), where A£' is the energy gap between the up and 
down states of the spin. Kb is Boltzman's coefficient and T is the tempera- 
ture of the thermal heat bath. Note that I\E = 2^B with B the magnetic 
field, and 7 the material dependant gyromagnetic constant which depends 
on the nucleus or particle, and is thus responsible for causing the differences 
in polarization biases [the ^^Carbon's nucleus, for instance, is 4 times more 
polarized than Hydrogen's nucleus; electron-spin is about 10^ times more 
polarized than Hydrogen's nucleus.] For high temperatures or small biases 



we approximate 



^ (6) 

2KbT ^ ' 



to leading order. A spin temperature at equilibrium is therefore T = 2KBta.nh-^{e) ^ 
which is T to leading order (for small e) . A single spin can be viewed 

as a single bit with "0" meaning spin-up and "1" meaning spin-down. Then, 
the single spin entropy (for a small polarization bias) is 

H = -Pi log2 Pt - n log2 n ~ 1 - (eV In 4) + O(e^) . (7) 

A spin temperature out of thermal equilibrium is still defined via the same 
formulas. Therefore, when a system is removed from thermal equilibrium, in- 
creasing the spins' polarization bias is equivalent to cooling the spins (without 
cooling the system) and to decreasing their entropy. What we have learned 
about single spins can be generalized to several spins on a single molecule. 
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One method for effective cooling of the spins is the reversible polarization 
compression (RFC). As is explained in the paper, RPC techniques essentially 
involve harnessing the powerful data compression tools which have been de- 
veloped over the past several decades and making use of them to cool spins, 
as much as is allowed by Shannon's bound on entropy manipulations in a 
closed system Let us present in more details the example of how Shan- 
non's bound limits RPC techniques. The total entropy of the n uncorrelated 
spins with equal biases satisfies H{n) = n(l — e^/ ln4). This entropy could 
be compressed into m > n(l — e^)/ln4 high entropy spins, leaving n — m 
extremely cold spins that have almost zero entropy. Due to the preservation 
of entropy, the number of extremely cold bits cannot exceed ne^/ln4. With 
e = 10~^ for instance, extremely long molecules whose length is of an order 
of magnitude of 10^° are required in order to cool a single spin close to zero 
temperature. If we use smaller molecules, with n <^ 10^'', and we compress 
the entropy into m = n — 1 fully-random spins, the entropy of the remaining 
single spin satisfies if (single spin) > n{l — / InA) — {n — 1) = 1 — ne^/ In 4. 
Thus we can at most reduce its entropy to 1 — [cg^g^j/ In 4] > 1 — ne^/ In 4, so 
that its polarization bias is improved to 



Unfortunately, manipulating many spins, say n > 100, is a very difficult task, 
and the gain of ^/n in polarization here is not nearly substantial enough to 
justify putting this technique into practice. It is interesting to note that 
improving the SNR by a factor of M can either be performed using n = 
spins via RPC, or using repetitions via repeated sampling over time, or 
using samples at once, via increasing the sample size times M^. 

The other effective cooling method is the special case of RPC called Po- 
larization transfer. Polarization transfer is limited as a cooling technique be- 
cause the increase of the polarization bias of the observed spins is bounded 
by the spin polarization bias of the highly polarized auxiliary spins. The 
polarization transfer method is commonly implemented by transfering the 
polarization among nuclear spins on the same molecule. This form of imple- 
mentation is regularly used in NMR spectroscopy, and it provides, say, an 
increase of about one order of magnitude. As a simple demostration, consider 
for example, the 3-bit molecule trichloroethylene (TCE) shown in Figure 3 
in the paper. The Hydrogen's nucleus is four times more polarized than each 
of the Carbons' nuclei, and can be used via polarization transfer to cool a 
single Carbon (say, the Carbon reffered to as C2), by a factor of four. 




(8) 
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A different form of polarization transfer involves removing entropy from 
the nuclear spins into electron spins or into other molecules. This technique 
is still at initial stages of its developement f^, but might become extremely 
important in the future as is explained in the Discussion of the paper and in 
Appendix lA.ll 

A.l Polarization transfer with electron spins 

Suppose that performing polarization transfer with spins other than the nu- 
clear spins on the same molecule (e.g., spins on other types of molecules, or 
electron spins on the same molecule) were to become feasible. It would clearly 
open interesting options potentially leading to a much more impressive effec- 
tive cooling than the regular polarization transfer onto neighboring nuclear 
spins. Polarization transfer with electron spins |lj could lead to three and 
maybe even four orders of magnitude of polarization increase. Unfortunately, 
severe technical difficulties with the implementation of this method, such as 
the need to master two very different electromagnetic frequencies within one 
machine, have so far prevented it from being adopted as a common practice. 
However, if such polarization transfer steps come into practice, and if the 
same machinery were to allow conventional NMR techniques as well, then 
algorithmic cooling could be applied with much better parameters. First, eg 
could be increased to around 0.01 or even 0.1. Second, the ratio -Rreiax-times 
could easily reach 10^ or maybe even 10^. With these numbers, scalable 
quantum computers of 20-50 bits might become feasible, as can be seen in 
the tables presented in Appendix |B| 

B Comparison with the original algorithmic 
cooling 

In order to compare our practicable algorithmic cooling algorithms with the 
BMRVV algorithmic cooling we need to consider longer molecules, and more 
time steps than considered in the text. This clearly makes the schemes unfea- 
sible with existing technology. However, if such parameters become feasible, 
for instance using the techniques described in Appendix lA.ll then medium- 
scale quantum computing devices of size 20-50 quantum bits can be built. 
It is important to mention that a comparison of the BMRVV algorithmic 
cooling with algorithims for RPC is provided in the supplementary material 
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Table 1: Sample results of the original algorithmic cooling [H] of 20 compu- 
tation bits (and with £ = 5). 



(^initial 


(^desired 


Jf 


^If 


T- 


eo < 1 




3 


140 


25x10^ 


16eo 


4 


180 


125x10^ 



(Appendix C) of 5J, to clearly show the advantage of the BMRVV algorith- 
mic cooling over RPC 

In the following, we use the structure in which there are n computation 
bits and a reset bit attached to each computation bit. For any desired final 
number of cooled bits m, the algorithm described in the text, namely PACl 
(Practicable Algorithmic Cooling 1), requires n = 2Jf + m computation bits, 
and the number of total time steps is mTj^. These numbers are obtained 
as trivial generalizations of Eq. 4 and 5 in the paper. We compare the 
performance of our new algorithmic cooling and the BMRVV algorithmic 
cooling |Sj in a case where the goal is to obtain 20 extremely cold bits. 

To obtain 20 cooled bits via Practicable Algorithmic Cooling 1 (PACl) 
we need to use n = 2Jf + 20 computation bits and T < 50 x 3'^-'"^ time 
steps This can be compared with j2j where n 40 jf [Eq.(9), with i = 5] 
and T < 400 x 5-^^+^ [Eq.(lO), with i = 5]. However, the required jf of jH] is 
different from the required Jf of the current work, because ej ~ 2ej_i in jS] 
while ej ~ (3/2)ej_i here (see Eq 2 in the paper). As a result, a smaller jf is 
required in the algorithm of 0, if a similar desired polarization bias is to be 
obtained. Tables and |21 present a fair comparison between the time-space 
requirements of the original and the new algorithms. As can be seen, the 
improvement obtained by our practicable algorithmic cooling (PACl) is very 
impressive, both in space and in time. 
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Table 2: Sample results of our improved algorithmic cooling of 20 computa- 
tion bits. 



(^initial 


(^desired 


Jf 






eo < 1 


7.6eo 


5 


30 


4040 


17.1eo 


7 


34 


36440 



C Quantum bits and quantum gates in NMR 
quantum computing 

Although here we use the language of classical bits, spins are actually quan- 
tum systems, and spin-half particles (two-level systems) are called quantum 
bits (qubits). A molecule with n spin-half nuclei can represent an n-qubit 
computing device. The quantum computing device in this case is actually 
an ensemble of many such molecules. In ensemble NMR quantum comput- 
ing [12] each computer is represented by a single molecule, such as the 
TCE molecule of Figure 3, and the qubits of the computer are represented by 
the nuclear spins. A macroscopic number of identical molecules is available in 
a bulk system, and these molecules act as though they are many computers 
performing the same computation in parallel. The collection of many such 
molecules is put in a constant magnetic field, so that a small majority of the 
spins (that represent a single qubit) are aligned with the direction of that 
field. To perform a desired computation, the same sequence of external pulses 
is applied to all the molecules/computers. Finally, a measurement of the state 
of a single qubit is performed by summing over all computers/molecules to 
read out the output on a particular qubit on all computers. In practicable 
algorithmic cooling we know in advance the polarization bias of the cooled 
bits after each step. In BMRVV algorithmic cooling we could also calculate 
the polarization bias, but it becomes a very cumbersome process. Thus, due 
to the complexity of the algorithm, the law of large numbers is used, and 
calculation of the polarization bias is done only after some particular steps 
called CUT in |]. 

For most purposes of this current theoretical work, the gates we consider 
are classical, and the spins can therefore be considered as classical bits. Only 
in a lab, the classical gates are implemented via the available quantum gates. 
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Then, the spins must be considered as qubits and the entire process is a simple 
quantum computing algorithm. However, unlike in other popular quantum 
algorithms, the use of quantum gates will not produce here any significant 
speed-up. 



D Details of new bias calculations 

as we have seen, the 3B-Comp can be built, for instance, via the following 
two gates: 

1. Use bit C as a control, and bit S as a target; apply a CNOT (ControUed- 
NOT) operation: C — > C, S — > S ® C, where ® denotes a logical 
eXlusive OR (XOR). 

2. Use bit B as a control, and bits A and C as targets; apply a variant of 
a C-SWAP operation: A ^ CB + AB, B ^ B, C ^ AB + CB; this 
means that A and C are swapped ii B — 0. 

The effect of a 3B-Comp step in case it is apphed onto three bits with equal 
bias e <^ 1 is as follows: 

• ii B — after the CNOT operation, the bias of C at that stage is 
(2e)/(l + e^) (this is the probabihty of S = C = given B ^ C, all 
these calculated prior to the CNOT operation; in other words, this is 
the probability of bit C being 0, given that bit B is after the CNOT). 
Due to the CSWAP this bias is transfered to bit A; 

• if i? = 1 after the CNOT operation, the bias of A is still e after the 
CSWAP. 



The new bias, enew) is the weighted average e 
X 2 



A ~ A 



26_ I 

2 ' 



1-62 



1+6 

e of the two possible biases which is enew — e 1 + , which 



4 

gives 

enew ~ 3e/2 (9) 

for small e. 

The following single operation can replace the CNOT+CSWAP, perform- 
ing a 3B-Comp, to cool bit A (the heated bits BC are not in the same state 
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as for the previous 3B-Comp, though): 



input : ABC output : ABC 

000 ^ 000 

001 ^ 001 

010 ^ 010 

011 ^ 100 

100 ^ Oil 

101 101 

110 110 

111 111 

Note that reset bits can be used for the computation as well. Thus the 
simplest algorithmic cooling can be obtained via 2 computing bits and one 
reset bit. 

Note also that the strong tools of data compression can easily be used 
to significantly improve the algorithms, for the price of dealing with more 
complicated gates. 



E A more space-efficient practicable algorith- 
mic cooling 

We consider here the case in which steps 1 and 2 (the 3B-Comp and PT steps) 
in the outline of algorithmic cooling are combined into one generalized step 
of RPC. Then, the logical gates are applied onto all bits in the system, that 
is computation and reset bits, to push the entropy into the reset bits. For 
comparison with the RPC and PACl (see Figure 1) we consider the cooling 
of a single bit. 

Once we use the reset bits for the compression steps as well, replacing 
the 3B-Comp+PT steps by a generalized RPC, wc can much improve the 
space complexity of the algorithm relative to PACl. We call this improved 
algorithm PAC2. 

Whenever we used 2n bits in PACl, namely n computation bits, each 
one having a reset bit as a neighbour, we can now use exactly n bits, namely 
n — 1 computation bits plus one reset bit. Let us explicitly show how this 
is done. Let eo be the polarization bias of the reset bit. In order to cool 
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a single bit, say bit A, to ei start with two computation bits, AB and one 
reset bit, C, perform PT(C B) followed by PT(i? ^ A to initiate bit 
A, RESET(C) (by waiting), then perform PT(C B) to initiate bit B, 
and another RESET(C). If the thermalization time of the compuation bits 
is sufficiently large, we now have three bits with polarization bias eo- Now 
perform 3B-Comp to increase the polarization of bit A. 

In order to cool one bit (say, bit A) to the second purification level (po- 
lariation bias £2), start with 4 computation bits (ABCD), and one reset bit 
E. Perform PT steps sequentially (with RESET(£') when needed) to initiate 
bits ABC, followed by a 3B-Comp on bits ABC that cools bit A to a bias 
€1. Then, perform PT steps sequentially (with RESET(£') when needed) to 
initiate bits BCD, followed by a 3B-Comp on bits BCD that cools bit -B to a 
bias ei. Next, perform PT steps sequentially (with RESET (ii^) when needed) 
to initiate bits CD, followed by RESET(i5), and then followed by a 3B-Comp 
on bits CDE that cools bit C to a bias ei. Finally, apply 3B-Comp to bits 
ABC to purify bit A to the second purification level €2 ~ (3/2)^eo. Clearly, 
the same strategy can be used to cool a single bit to ej^ ^ (3/2)-' ■^'eo using 
2Jp computation bits and a single reset bit. The total space complexity is 
therefore 

2Jf + 1 (10) 

Note that (3/2)"^ = 81/16 which is slightly larger than 5, leading immediately 
to the numbers presented in Figure 1. Cooling by a factor of 25 requires 17 
spins here, while it requires 625 spins in RPC, or it requires repeating the 
experiment 625 times if no cooling is used. 

We should mention, though, that the timing considerations for PAC2 
are more demanding than the timing consideartions for PAC2 due to the 
need of many more SWAP gates. All the steps of the algorithm must be 
done within the relaxation time of the computation bits, and even more 
demanding, within the dephasing time (known as T2) of the computation 
bits. 
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